Orbital and spin relaxation in single and coupled quantum dots 
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Phonon-induced orbital and spin relaxation rates of single electron states in lateral single and 
double quantum dots are obtained numerically for realistic materials parameters. The rates are 
calculated as a function of magnetic field and interdot coupling, at various field and quantum dot 
orientations, ft is found that orbital relaxation is due to deformation potential phonons at low 
magnetic fields, while piezoelectric phonons dominate the relaxation at high fields. Spin relaxation, 
which is dominated by piezoelectric phonons, in single quantum dots is highly anisotropic due to 
the interplay of the Bychkov-Rashba and Dresselhaus spin-orbit couplings. Orbital relaxation in 
double dots varies strongly with the interdot coupling due to the cyclotron effects on the tunneling 
energy. Spin relaxation in double dots has an additional anisotropy due to anisotropic spin hot 
spots which otherwise cause giant enhancement of the rate at useful magnetic fields and interdot 
couplings. Conditions for the absence of the spin hot spots in in-plane magnetic fields (easy passages) 
and perpendicular magnetic fields (weak passages) are formulated analytically for different growth 
directions of the underlying heterostructure. It is shown that easy passages disappear (spin hot 
spots reappear) if the double dot system loses symmetry by an xy-like perturbation. 

PACS numbers: 72.25.Rb, 73. 21. La, 71.70.Ej, 03.67.Lx 



C^ 



X3 

o 



> 

m 

en 

o 
o 



I 

o 



X 



I. INTRODUCTION 

The spin degree of freedom in solid state systems has 
a much longer memory than orbital degrees. This fact is 
exploited in spin electronicsi as well as in spin quantum 
computing, most notably in quantum dots- in which elec- 
tron spin provides a qubit for controlled single and two- 
qubit operations2i^. The performance of the spin qubits 
is ultimately limited by spin relaxation and decoherence. 
At present we believe that general principles and mecha- 
nisms of spin relaxation and decoherence in quantum dots 
arc known, while it remains to develop particular models, 
understand realizations of the mechanisms as well as to 
perform realistic calculations, in special cases of interest. 

There appear two principal mechanisms of spin relax- 
ation in quantum dots. At low magnetic fields (say, 
tens of gauss), the relaxation proceeds via hyperfine 
coupling of the electron spin with the lattice or im- 
purity nucleii^*^ On the other hand, at higher fields 
(Teslas) the relaxation is due to phonon-induced spin- 
flip transitiong7i8i9d"iiid2d3d4ii5,i6,i7,i8,i9_ ^^^ese are al- 
lowed due to the presence of spin-orbit coupling. Vari- 
ants of the phonon-induced spin relaxation has been 
proposed, such as ripple coupling, ^"'^^ important in 
very small quantum dots (10 nm), or fluctuations in 
spin-orbit parameters, important when underlying het- 
erostructure inhomogeneitiesS^ are present. A possi- 
ble direct spin-phonon couplingSiSiS^, due to spin-orbit 
modulated electron-phonon coupling, have been found 
of less importance. Phonons can also change the spin 
precession rate and cause spin decoherence.-^ As was 
shown in Ref. .25 phonon-induced spin relaxation and 
decoherence proceed on similar time scales. Another 
source of the decoherence is the fluctuation of the gate 
potential.^ The experimental results on spin relaxation 
in single^i22i224^Si2LS and double dots^SiM, as well as 
on orbital relaxationSi2&, support the above theoretical 



picture. 

Here we present a systematic and comprehensive in- 
vestigation of phonon-induced orbital and spin relax- 
ation in lateral single and double quantum dots, defined 
in a GaAs heterostructure. We consider the most rel- 
evant electron-phonon couplings — the deformation po- 
tential and piezoelectric acoustic phonons, and realistic 
spin-orbit couplings — the Bychkov-Rashba and the lin- 
ear and cubic Dresselhaus ones. We calculate the relax- 
ation rates in the presence of in-plane and perpendicu- 
lar magnetic fields. Our results are numerically exact 
within the Fermi's Golden rule approximating the tran- 
sition rate. We have already reported on new anisotropy 
effects of spin relaxation in double dots, in Ref. [S^. The 
anisotropy arises due to anisotropic spin hot spots, the 
parameter (magnetic field and interdot coupling) regions 
in which a spectral crossing between a spin up and a 
spin down state is lifted (producing an anti-crossing) by 
spin-orbit coupling22i2S. At these points the spin and or- 
bital relaxation rates are equal. In single quantum dots 
spin hot spots were found in Ref. 40 . while in vertical 
few-electron quantum dots in Ref. [43. In lateral double 
dots spin hot spots appear at useful magnetic fields (1 
T) and interdot couplings (0.1 meV), due to the cross- 
ing of the lowest orbital antisymmetric (with respect to 
the quantum dot axis) state and the Zeeman split sym- 
metric state of the opposite spin42i This occurs when the 
tunneling energy equals the Zeeman energy. Manipula- 
tion of interdot coupling in the presence of a magnetic 
field thus in general results in a short spin lifetime. For- 
tunately, we have found^^ that the spin hot spots are 
absent for certain arrangements of the double dots' axis 
and the orientation of the in-plane magnetic field. In 
particular, if the dots arc oriented along a diagonal [on 
a (001) heterostructure plane], and the magnetic field is 
perpendicular, the spin hot spots are absent (due to sym- 
metry reasons) for any values of spin-orbit parameters. 



We have proposed such a geometry, which corresponds to 
what we cah "easy passage"—, for quantum information 
processing. 

The particular resuhs of our Letter, Ref. l23j are not 
repeated here. Instead we focus on providing a unified 
description, both analytical and numerical, of orbital and 
spin relaxation rates. We give analytical formulas de- 
scribing the trends, with respect to magnetic fields and 
confinement energies of the dots, of the rates. We present 
the numerically calculated orbital relaxation rates and 
demonstrate that they are due to the deformation poten- 
tial phonons at low magnetic fields and due to piezoelec- 
tric phonons at high fields (at zero magnetic field orbital 
relaxation in a biased double dot was studied in Ref. Il3 . 
using a two-level model). 

As for spin relaxation, we demonstrate here the differ- 
ent origin of spin hot spots in single and double quantum 
dots. While in single dots spin hot spots appear due to 
the Bychkov-Rashba coupling'*^, in double dots both the 
Bychkov-Rashba and Dresselhaus couplings contribute. 
The reasons is the different symmetry of the underly- 
ing states in single and double dots^^. Furthermore, we 
classify here the conditions for the absence (or narrow- 
ing) of spin hot spots in double dots defined in quan- 
tum wells grown in different crystallographic directions, 
in which the Dresselhaus spin-orbit coupling takes on dif- 
ferent functional forms. We also explore the orbital ef- 
fects of a perpendicular magnetic field component — the 
main effect is the absence of easy passages; only a nar- 
row "weak passages" appear instead with inhibited but 
finite spin hot spots. Finally, we show that easy passages 
are also absent in general asymmetric double dots, im- 
plying stringent symmetry requirements on coupled dots 
for spin information processing. 

The paper is organized as follows. In Sec. II we de- 
scribe our model of single and double quantum dots, de- 
rive relevant perturbations responsible for the spin relax- 
ation, and write useful expressions for orbital and spin 
relaxation rates. In Sec. Ill we describe the orbital and 
spin relaxation relaxation in single dots for the case of 
in-plane and perpendicular magnetic fields. Section IV 
gives a similar description for double dots. Finally, in 
Sec. V we give conclusions. 



II. MODEL 

A. Hamiltonian 

We study a two-dimensional electron gas confined in a 
(001) plane, spanned along x and y directions, of a zinc- 
blende semiconductor heterostructure. An additional 
confinement into lateral quantum dots is given by top 
gates. The single electron Hamiltonian, in the presence 
of magnetic field and spin-orbit coupling, is 



Here T = P^/2to is the kinetic energy with the effective 
electron mass m and kinematic momentum P — p-|-eA = 
—ih\I + eA; e is the proton charge and A is the vector 
potential of the magnetic field B. If the magnetic field is 
perpendicular to the plane, B^ = (0,0,5^), we choose 
the vector potential as A^ = (i?^/2)(— y,2;,0). If the 
field is in the plane, lying under the angle 7 relative to 
X, B|| = _B||(cos7,sin7,0). The orbital effects of the in- 
plane field are inhibited — only the Zeeman interaction is 
taken into account in this case. The position vector is 
denoted as R = (a;, y, z) = (r, z). We will also find useful 
to introduce the kinematic angular momentum L = R x 
P. 

The quantum dot geometry is defined by the confining 
potential 

V{r) = -mujlmm{{r - df, (r + d)^}. (2) 

The distance between the minima of the potential is 2d, 
which is further called the interdot distance. The an- 
gle between d and x is denoted as 6. If d=0, the po- 
tential is parabolic, V — [1 / 2)mu}'^r'^ , representing the 
single dot case with the confinement energy Eq = fkvo 
and the confinement length Iq — {h/mujo)^''^, setting 
the energy and length scale, respectively. Both the con- 
fining potential and the vector potential of the perpen- 
dicular magnetic field define the effective length Ib = 
lo{l + Ble^ll/Ah^)-^/^. 

The Zeeman energy is given by Hz — (.g/2)/iB(T.B, 
where g is the conduction band g factor, ^b is the Bohr 
magneton, and cr are the Pauli matrices. To shorten 
the notation, we will use a renormalized magneton /i = 

{g/2)fiB- 

Spin-orbit coupling gives three important terms in con- 
fined systemsi The Bychkov-Rashba HamiltonianjiSili 



HbR = OIBB. {c^xPy - CTyPx) /K 



(3) 



appears if the confinement is not symmetric in the growth 
direction (here z). The strength cilbr of the interac- 
tion can be tuned by modulating the asymmetry by top 
gates. Due to the lack of spatial inversion symmetry in 
zinc-blende semiconductors, the spin-orbit interaction of 
conduction electrons takes the form of the Dresselhaus 
Hamiltonian^ which gives two terms, one linear and one 
cubic in momentumii& 



Hd = IciP^/n^) {-(TxPx + (JyPy) /K 



(4) 



H^T + V+Hz + H, 
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Hd3 = M2){a,PxP^^ayPyP^)/h^ + H.c., (5) 

where jc is a material parameter. The angular brack- 
ets in Hb denote quantum averaging in the z direction- 
the magnitude of Hb depends on the z confinement 
strength. We express the strengths of the linear spin- 
orbit interactions in length units by Ibr — fi^/2maBR. 
and Ib = h'^/2mjc{P^). 

In our numerical computations we use bulk GaAs ma- 
terials parameters: m = 0.067 Trie, g = —0.44, and 



7c = 27.5 eVA^. For the coupling of the linear Dres- 
selhaus terms we choose ^c^Pl) j'f? ~ 4.5 meVA, corre- 
sponding to the 11 nm thick ground state of the trian- 
gular confining potential4i To agree with experimental 
data (see Ref. i32i) we choose for (Xbr a value of 3.3 
meVA, which is in line of experimental observations^SiiS 
and corresponds to the carrier density of 5x10^^^ cm"^ 
in Ref. !50|. These values correspond to length units of 
Ibr — 1-8 /xm, and lu — 1.3 /xm. 

For a confinement length of 32 nm (used in a re- 
cent experiments^) and a perpendicular magnetic field 
of 1 T, one gets the following typical magnitudes for the 
strengths of the contributions to the Hamiltonian given 
by Q: 1-1 nieV for the confinement energy £'o, 13 /ieV 
for the Zeeman splitting, and 14, 10, and 0.8 /j,eV for 
the linear Dresselhaus, Bychkov-Rashba, and the cubic 
Dresselhaus terms, respectively. The spin-orbit interac- 
tions are small perturbations, with strengths comparable 
to the Zeeman splitting. This leads to the many orders 
of magnitude difference between the orbital and spin re- 
laxation rates. 

We numerically diagonalize the full Hamiltonian (Q 
(see Ref. |43 for further details) and compute the orbital 
and spin relaxation rates using Fermi's Golden Rule. We 
also present analytical calculations for various limiting 
cases. 



B. Perturbative eigenfunctions 

Our numerical results can be qualitatively understood 
from considering the lowest order of the perturbative so- 
lution of the Hamiltonian Q. We assume that spin- 
orbit couplings are small perturbations and that one can 
solve the Schrodinger equation for Hamiltonian Hq = 
T + y+ Hz- First, we transform^- the Hamiltonian to 
remove the linear spin-orbit terms 
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Ho + Hi, 
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(7) 



The transformation is defined by operator iJ"^ = -ff^'^ + 

H°P = {i/2lBR){yax - xoy) - {i/2lD)ixax - yay), (8) 

Keeping only terms up to the second order in the linear 
spin-orbit and Zeeman couplings and the lowest order 
term in the cubic Dresselhaus coupling, the new terms of 
the transformed Hamiltonian are 



H^^^ = [HD.H°ji']/2 = -{te/Amll){l~L,<y, 



(9) 



r(2) 



H);k = [Hbr. i?S,]/2 = -{h'/Aml'BR)(l + i.a,ij;0) 

+ {nB\\/lD)(Tz{xsu\-f + ycos-f), (11) 

Hbr ^ [Hz, H°^r\ = {lJ.B±/lBR){ycry + xa^) - 

-{^iB\\/Ibr)(Jz{xcos^ + y&m^). (12) 



(2) (2) 

While Hj^ and Hg^^ are transformed linear spin-orbit 
terms in the absence of magnetic field, the terms H^ 
and -ff^^ describe the mixing of the spin-orbit and Zee- 
man interaction in the unitary transformation given by 
H°^] these terms are essential for understanding spin re- 
laxation anisotropy. 

We denote the eigenfunctions and eigenenergies of i?o 
as ^ and e. We use the standard perturbation theory 
for non-degenerate states and then remove the unitary 
transformation to get the approximate eigenfunctions, ^, 
of the original Hamiltonian |^ : 



*. 






(13) 
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For degenerate states, which normally lead to spin hot 
spots with strongly enhanced spin relaxation, we use 
Lowdin's perturbation theory. If two eigenstates, ^^ and 
^I'j, of Hq are degenerate, the corresponding perturbed 
states are 



*,: = e-^"' 



Pa ^ i 



Pij w j 



E 



f-i - e/c 



*A 



(14) 

The coefficients (3 are the solutions of the appropriate 
secular equation. If e^ — €j ^ {Hi)ij, then Pa ^ 1 
and (3ij <C 1 - Eq. \i'A\ is recovered. The other case, 
when Ci — ej < (Hi)ij, describes anti-crossing — spin hot 
spotsi^2i2S In the limiting case, when e^ = ej, we get 

A, = l/V2 = /3y. 

The analytical solution of Hq for the single dot case is 
known. The eigenstates are the Fock-DarwinSiiS states, 
^n,i,cr, where n is the principal quantum number, I is 
the orbital quantum number, and a is the spin. For the 
double dot case the analytical solution is not known, but 
for our double dot potential the eigenfunctions can be 
approximated by a properly symmetrized linear combi- 
nation of Fock-Darwin functions centered at the potential 
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C. Phonon-induced orbital and spin relaxation 
rates 

By orbital relaxation we mean the transition from the 
first excited orbital state to all lower lying states. By 
spin relaxation we mean the transition from the upper 
Zeeman split orbital ground state to all lower lying states 
(except at high magnetic fields, there is only one lower 
Zeeman split orbital ground state). The spin of a state ^ 
is quantized in the direction of the magnetic field. How- 
ever, due to the spin orbit interactions, the perturbed 
states ^ have no common spin quantization axes. We 
call a state to be spin up (down) if the mean value of 
the spin in the direction of the magnetic field is positive 
(negative). Since the spin-orbit interactions are a small 



perturbation, these mean values are close to ±1/2, except 
at anti-crossings. 

Given the initial and final states for the transition we 
compute the rates by Fermi's Golden rule, where the per- 
turbation is the electron-phonon interaction. The rele- 
vant terms for our GaAs system comprise deformation 
(df) and piezoelectric acoustic (pz) phonons, described 
by Hamiltonian terms^^ 



H 



df 



(TeY, 



K 



' hK 

2pVci 



{b^.i 






„iK.R 



(15) 



Here the three-dimensional phonon wave vector is de- 
noted by K = {kx,ky,kz) — (k, /c^), and A ~ l,tl, or 
t2 is the phonon's polarization (one longitudinal and two 
transversal); p is the material density (5.3 x lO'^kg m^'^, 
for GaAs), V is the volume of the crystal, c\ is the 
sound velocity, (c; = 5.3 x 10^ m/s, ct — 2.5 x lO'^ m/s), 
^K A'^K,A are the creation and annihilation phonon op- 
erators, Ue is the deformation potential (7.0 eV), and 
ehi4 is the piezoelectric constant (1.4 x 10^ eV/m). Fi- 
nally, the geometrical factors M\ are equal to 2{kxkye^ + 
kzkxSy +kykze^) / K'^ , where e^ arc unit polarization vec- 
tors. 

Consider first the deformation potential, Eq. 1151 in 
which only longitudinal (A = I) phonons take part. Using 
Fermi's Golden rule, a relaxation (orbital or spin) rate 
can be written as 



-^dt 



n{E) + l]j^,E' / d2k|F(k)n/(fci)|Vfci(17) 



(18) 



[n{E) + l]jdiiEyiB)xdfm. 



Here E is the energy difference between the initial and 
final states, n{E) = [exjp{E/kBT) — 1]^^ is the oc- 
cupation number of the phonon state with energy E 
at temperature T (further we use zero temperature), 
7dt = a-^/STT^pcfh'^ is the strength of the deformation 
electron-phonon interaction [8.3 x 10^*^ s~^ nm /(meV)^], 

F(k) = /d^rfjC*'^'"!'/ is the xy-overlap, and /(fc^) = 
/ dz iple^'^'^^^po is the z-overlap, contribution of which can 
be neglected, /(fcz) w 1, if the energy difference E is 
much smaller than the excitation energy in the z confine- 
ment potential. The z component of the wave vector is 
given by k^ = \/£'x/Ib ~ ^^' where the dimensionless pa- 
rameter £x = Els/hcx is the ratio of the effective length, 
Zb, and the the wavelength of the emitted phonon. Fi- 
nally, Xdi{£i) is an integral of the xy-overlap i^(k). Since 
the typical linear dimension of a wave function is the ef- 
fective length Ib, we express it as 



Xdf(fi) 



fc,>o ^j£i-kH 



(19) 



We compute the relaxation rate F'*^ numerically using 
formula (|17|) . However, we can gain physical insight in 
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TABLE I: The relaxation rates and the relative strength of 
the contributions due to deformation (A = I), piezoelectric 
longitudinal (A = I), and piezoelectric transversal (A = t) 
phonons. The two limiting cases are defined by the ratio, £x, 
of the wavelength of the emitted phonon of polarization A, 
and the effective length Ib- The initial and final states are 
encoded into the coefficient A, which needs to be evaluated 
for specific cases. 



two important limits. First, if the wavelength of relevant 
phonons is smaller than the size of the dots, £i ^ 1, 
the square root can be taken out from the integral and 
Xdf ^ 1/^/- Physically, this means that the energy to be 
absorbed by the phonon is large and almost whole is in 
the z component of the phonon wave vector (phonon is 
emitted almost perpendicularly to the xy-plane). Sec- 
ond, in the opposite limit oi £i <C 1, the integration 
is only in the vicinity of point k = 0. Because of 
the orthogonality of the eigenfunctions the overlap in- 
tegral vanishes, i^(k — *■ 0) ^ 0, and the lowest order 
gives |Fp ~ (Mb)^- This leads to the dependence of 

Xdf{£i)^£l 

Analogous expression holds for the piezoelectric in- 
teraction which contains contributions from longitudinal 
and transverse phonons. The relaxation rate can be writ- 
ten as 



rP''^[n{E) + l]jp,{l/lB)Xpz{£), 



(20) 



with 7pz = {ehii)'^/8TT^pcfh = 4 x 10^" s'^ (note 
the different unit from 7df) nm and Xpz(^) = 
Exi^f/cl)!dH\^lB)\Mx\'\F{k)\yk^lB. The geometri- 
cal factors, Mx, have no influence on the limiting expres- 
sions for Xpz in the limit £ <C 1, where Xpz(^) ^ £^- If 
£ ^ 1, the fact that Mx contains factors (k^/K)'^ and 
kx/K leads to limits £~^ and £~^ for the longitudinal 
and transverse phonons, respectively. Table ^ summa- 
rizes the limiting expressions. 

In addition to the deformation and piezoelectric 
phonons, there are additional electron-phonon spin de- 
pendent interactions which can lead to spin relaxation. 
However, a direct spin-phonon coupling (spin-orbit mod- 
ulated electron-phonon interaction^'^) is believed to give 
a negligible contribution. In very small (say, 10-20 nm, 
which is not our case) quantum dots spin relaxation due 
to the so called ripple mechanismSS can be as important 
as the spin-orbit mechanism and should be considered. 
Finally, at low magnetic fields the relaxation is believed 
to be dominated by the hyperfine interaction between the 
electron and nuclei of the host materiali^ 



III. SINGLE DOTS 

In the single dot case we identify the unperturbed lower 
and upper Zeeman split orbital ground, and excited or- 
bital states as \E'o,o,Tj ^o,o,ii and \l/o,-i,t, respectivelyi^ 
The negative value of the g factor energetically favors 
spin up rather that spin down states. Having opposite 
spin, the perturbed ground and spin states will have a 
nonzero overlap due to those perturbations in the trans- 
formed Hamiltonian ® which do not commute with the 
Zeeman term. Therefore the xy-overlap F(\i) will be 
proportional to the strengths of the corresponding per- 
turbations. In the case of the spin relaxation, the coef- 
ficient A in Tab. ^will be approximately equal to these 
strengths divided by a typical energy difference between 
the corresponding coupled states, as can be seen from 
Eq. H13|l . On the other hand, since the excited orbital 
and ground states have the same spin, the coefficient A 
for the orbital relaxation is of order 1. 

This consideration leads to the following approxima- 
tions which we use when estimating the rate analytically. 
For orbital relaxation 



Hi wO. 



(21) 



For spin relaxation, in analytical calculations we neglect 
the cubic Dresselhaus term. If the magnetic field is in- 
plane, Eqs. lO and l|10() do not couple the ground or 
the spin state, which have zero orbital momenta, with 
any other state. If the magnetic field is perpendicular, 
Eqs. © and (|10|l commute with the Zeeman term, again 
giving no contribution to the spin relaxation. For the 
spin relaxation we therefore approximate 



Hi 



Hd + Hbr- 



(22) 



A. In-plane magnetic field 



We have earlier^^ compared our calculation of the 
spin relaxation in an in-plane magnetic field with the 
experiment S We have found that the experiment can be 
explained with a reasonable set of spin-orbit parameters 
which we use also in the present article. Using Eqs. (|13|l . 
(|20|l , and (|22() we get for the dominant contribution to the 
spin relaxation due to piezoelectric transversal phonons 
in the low magnetic field limit 
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where 



ll, + llR-2sin{2-f)lDlBR 



-so 



''D'-BR 



(23) 



(24) 



describes the effective (anisotropic) spin-orbit length. 
The angular dependence of the spin relaxation rate, ex- 
pressing the C2v symmetry of the heterostructure, allows 



to find the ratio of the Dresselhaus and Bychkov-Rashba 
couplings: 



Ta:irL{lD/lBR,lBR/lD} = 2/(^r^+ 1) - 1, 



(25) 



where r^ is the ratio of the rates at 7 = 45° and 7 = 
135°. A possible measured angular dependence with the 
minimum at 7 = 45° would be a convincing indication 
that the admixture due to spin-orbit is the mechanism 
of the relaxation. A more general angular dependence, 
allo wing for out-of-plane magnetic fields, was derived in 

Ref. Ila 

The reason for the angular dependence of F^^ follows 
from Eq. (|22|) . which for an in- plane field is 



Hi = -^jLB\\a:,[x{ 
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1 
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sm7 sm7 

-)+y(- 
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BR 



cos 7 



)]• (26) 



Due to the selection rules for the Fock-Darwin states, x 
and y do not mix in coupling of the states. The coeffi- 
cient A^ is then proportional to the sum of the squared 
couplings from Eq. H2()|) . at x and y. Taking Eq as a typi- 
cal energy difference E of the coupled states and using Is 
for a natural length unit, we get A^ w \fiB\\lB/Eo\'^CgQ. 
Noting that Ib = lo for in-plane field and using the low 
energy limit for F^^ from Tab. HJ one recovers Eq. H23|) 
up to a numerical factor. The numerical result has been 
presented in Ref. |33 (Fig. 1) and is not repeated here. 



B. Perpendicular magnetic field 

1. Orbital relaxation rates 

In the case of a perpendicular magnetic field, the nu- 
merically calculated orbital and spin relaxation rates in a 
single dot are shown in Fig. ^ The orbital relaxation rate 
is of the order of 10^ s^^. The spin-orbit contributions to 
the rate (not shown in the figure) are of the order of 10^ 
s^^ for the linear spin orbit terms and 10^ s~^ for the cu- 
bic Dresselhaus term, validating the approximation Eq. 
(|21|l . The energy difference of the orbital and the ground 
state is £' = h?/2mlg — (he/2m)B±. At low magnetic 
fields the high £ limit applies and the deformation poten- 
tial dominates the orbital relaxation rate. The results are 
listed in Tab. ^ The values at zero magnetic field, up to 
a numerical factor, follow from Tab. Q if one uses A = 1 
and the low magnetic field limits, where E « ft^/mZ§, 
and Ib ~ Iq- The dependence of the rates on the energy 
difference of the states, shown in Tab. U is enough to 
understand the different dependence of the deformation 
and piezoelectric contributions to the orbital relaxation 
rate at low magnetic fields shown in Figs, ^andj^l The 
deformation contribution drops with increasing both the 
magnetic field and confinement lengths, while the piezo- 
electric contribution increases with increase of these two 
parameters. 

For fields lower that 1 T the dominant deformation 
contribution manifest itself on Fig. |21 At magnetic fields 




T 6 

magnetic field B [T] 



10 



FIG. 1: (Color online) Orbital and spin (labels D,BR, and 
1)3 denote which spin-orbit interaction is present) relaxation 
rates in a single quantum dot, for the piezoelectric (solid) 
and deformation potential (dashed) phonons. The confining 
length is 32 nm. Anti-crossing of the unperturbed spin and 
orbital state occurs at B±_ = 5.2 T. 
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TABLE II: Approximate orbital and spin (due to Dresselhaus 
coupling) relaxation rates in a single quantum dot at low and 
high magnetic fields in lowest order of the non-degenerate 
perturbation theory. In the last column we state the maximal 
or minimal magnetic field by requiring that at Zo = 32 nm the 
presented approximation does not differ from the numerical 
value by more than a factor of 2. 



higher tlian 1 T the piezoelectric contribution dominates. 
Up to about 4 T we are still in the regime £ 3> 1 and the 
rate grows with increasing magnetic field and increasing 
confinement length. Since the energy difference E drops 
with increasing magnetic field, for magnetic fields ^ 6 T 
we get into the limit £ -^ 1. The corresponding orbital 
relaxation rates in^Jthen follows from Tab. fusing A — 1 
and high magnetic field limits, where E w ft' /emBl'^ and 
Z^ « 2h/eB. This leads to a much stronger drop of the 
deformation contribution to the rate with the increase of 
both magnetic field and the confinement length, than is 
the drop of the piezoelectric contribution. 
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FIG. 2: (Color online) Orbital relaxation rate (the sum of the 
deformation and piezoelectric contribution) in a single quan- 
tum dot 'as a function of magnetic field and the confinement 
length Iq I the confinement (excitation) energy i?o. The rate 
is given on the logarithmic scale in the units of s~^. The solid 
lines represent equal relaxation rates (equirelaxation lines). 
The granular structure in the figure is an artifact of the lim- 
ited data resolution. 



Finally, we explain the influence of the anti-crossing 
on the orbital relaxation rate, seen in Fig. ^ The anti- 
crossing contributes by an overall factor of |/3mP, see Eq. 
(|14|l . which multiplies the orbital relaxation rates listed 
in Tab. ^] Solving the appropriate secular equation, we 



gel, 



lA. 



E 



2^W^\C\ 



2' 



(27) 



where E = t? j2ml\ — {he/2m)B±^ — 2fj,B± is the en- 
ergy difference between the crossing states, and C = 
"mTl — (^ ~ B±el^/2h) is the strength of the coupling 
between these states due to the Bychkov-Rashba term. 
Away from anti-crossing Pa w 1, while directly at the 
anti-crossing the rate is reduced by a factor of 2. The 
anti-crossing region for the orbital relaxation is rather 
narrow {'^ 0.1 T) and manifests itself as a narrow line of 
the suppression of the rate in Fig. [5] 



2. Spin relaxation rates 

For spin relaxation the relevant energy difference is the 
Zeeman splitting, E w 2\fiB\. Therefore the low energy 
limit, £ <C 1, applies up to rather high magnetic fields 
(~10 T). Piezoelectric transversal phonons dominate the 
rate. The linear spin-orbit terms dominate over the cu- 
bic Dresselhaus term, although the difference becomes 
smaller for higher magnetic fields. We use an example 
of the linear Dresselhaus term for analytical expressions. 



Using Eq. (|22() and the limits of low and high magnetic 
fields we present the analytical spin relaxation rates in 
Tab. |n] (these results were also derived in Refs. Ho and 
Q). These formulas approximately follow from Tab. ^ 
using A — \ij.B±\Ib/IdSE and E = \iiB±\, while noting 
that SE = Eq for low and SE = \i.J.B±\ for high magnetic 
fields. The trends described by the Dresselhaus contri- 
bution can be seen in Fig. ^ The spin relaxation rate 
grows much steeper with increasing magnetic field at low 
B± (fifth power) than at high B± (first power). Interest- 
ingly, at high magnetic fields the rate does not depend 
on the confining length. 

Away from anti-crossing analogous formulas, up to a 
numerical factor, as those listed in Tab. ^ hold for the 
contribution to the spin relaxation due to the Bychkov- 
Rashba term after the substitution lo — *■ Ibr- In this 
case the contribution to the overlap between the spin 
and ground states due to the term f3ij in Eq. !(T^ is 
negligible. However, comparing the analytical formulas 
from Tab. ^]with the numerical calculation in Fig. |31 we 
find discrepancy, except at low magnetic fields. This is 
because, as can be seen also in Fig. ^ the rate is actually 
dominated by a spin hot spot (anti-crossing). The anti- 
crossing occurs for single dots only when the Bychkov- 
Rashba term is present, since the Dresselhaus terms do 
not couple the unperturbed orbital states42i^ In this 
case we can neglect all terms but that one containing 
Pij in Eq. (|14|) and for the spin relaxation rate due to 



the anti-crossing one gets r(spin, acr) 
The secular equation gives 



lA.f 



E 



2^WTW 



l^yl^r(orbital). 



(28) 



where the variables are those defined under Eq. (|27|l . 
Thus, the anti-crossing effectively mixes what we usually 
call spin and orbital rates. The spin relaxation rate has a 
sharp peak at the anti-crossing. With increasing the "dis- 
tance" from the anti-crossing the rate drops, mirroring 
the drop of the coefficient |/3ijp. Only far enough from 
the anti-crossing the term f3ij is negligible in Eq. (|14|l and 
the rate is described by expressions analogous to those 
from Tab. ^ In Fig. ^the Bychkov-Rashba contribution 
to the spin relaxation rate is dominated by the f3ij term 
unless the magnetic field is smaller that 2 T. Similarly in 
Fig. [3 for fields higher than 2 T the total spin relaxation 
rate is dominated by the anti-crossing contribution due 
to Bychkov-Rashba term. Consequently, the influence of 
the anti-crossing is substantial in a much larger region 
(several Tesla) than in the case of the orbital relaxation. 

In Ref. H^ spin relaxation rates due to the deforma- 
tion potential were computed in the lowest order of the 
perturbation theory and an analogous figure to our Fig. 
n was presented. Our results for both orbital and spin 
relaxation rates are in a quantitative agreement. 
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FIG. 3: (Color online) Spin relaxation rate in a single quan- 
tum dot as a function of magnetic field and the confinement 
length lo I the confinement energy Eq. The rate is given 
on the logarithmic scale in the units of s~^. The solid lines 
represent equirelaxation lines. 



IV. DOUBLE DOTS 



In our double dot potential the ground (excited orbital) 
state can be approximated as a symmetric (antisymmet- 
ric) combination of two Fock-Darwin functions, ^o,o,Tj 
placed at the two potential minima. In Ref. |42 we have 
studied the energy spectrum and classified the symme- 
tries of the states of a double dot with a potential given 
by Eq. (0). What we call here ground, spin and orbital 
state is denoted there as F^, F^, and F]j, respectively. 
The upper index indicates spin and the lower index in- 
dicates the symmetry of a particular state with respect 
to spatial inversion. The energy difference between the 
ground and excited orbital state, 26 Et, is strongly in- 
fluenced by the ratio of the interdot distance and the 
effective length^iS, V = d/ls , 



25Et^ 



h^ 2V{l-9^){l+V^EriciV)-e-'^'} 



(29) 

where a dimensionless parameter 9 = B±el%/2h. The 
tunneling energy 6Et gives the frequency of single- 
electron coherent oscillations between the left and right 
dots. The approximation of Eq. H22|l for spin relaxation 
is correct also here, since Eqs. (O and H1U|) do not cou- 
ple any two of the ground, spin, and orbital states due 
to a deflnite symmetry of the Lz operator42i There is a 
coupling through higher excited states with appropriate 
symmetry, but, as we learn from numerics, this is negli- 
gible. 
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FIG. 4: (Color online) Spin relaxation rate in a double quan- 
tum dot as a function of in-plane magnetic field for 7 = 0° 
and the interdot distance d / tunneling energy 5Et, for a con- 
finement length 32 nm. The relaxation rate is given on the 
logarithmic scale in the units of s~^. The double dot is ori- 
ented along [100] (5 = 0°). 



A. In-plane magnetic field 

The spin relaxation rate as a function of in-plane mag- 
netic field and the interdot distance is plotted in Fig. 
0] The rates for small interdot distances are similar to 
the single dot case, where the rate grows with increasing 
magnetic field; for low magnetic fields more steeply than 
for large. The order of magnitude of the rate is given by 
Eq. (ESJl, being about 10^ s'^ at 1 T and 10^ s'^ at 10 T. 
At large interdot distances the rate is strongly influenced 
by the presence of a anti-crossing (spin hot spot), which 
occurs when the Zeeman and twice the tunneling energies 
are equal.-^ If the tunneling energy is changed from zero 
to a value of order of the single dot excitation energy, re- 
gardless of the magnetic field strength, one always passes 
through a spin hot spot region, where the spin relaxation 
is very fast. Fortunately there exist specific orientations 
of the double dot system and the magnetic field, where 
this anti-crossing does not occur. We call such a config- 
uration "easy passage. "-fii 

To understand the angular dependence of the rate pre- 
sented in Ref. |33 and find conditions for an easy pas- 
sage we transform the Hamiltonian J^J) with Hi given by 
Eq. (|22|l into coordinates in which the new x axis lies 
along the dot's axis d. Since there arc no orbital effects 
in in-plane magnetic fields, in these new coordinates the 
unperturbed solutions of the Hamiltonian Hq have a def- 
inite symmetry under inversions about x - the ground 
and spin states are symmetric, while the orbital state is 
antisymmetric. The transformed Hi of Eq. H22I) , is 



In the single dot case the coefficient A^ in Tab. Q] is 
proportional to the sum of the squared couplings in Eq. 
%M)\ at X and y. However, in the double dot case, x 
and y can couple states differently. For large interdot 
distances the most important influence on the spin re- 
laxation comes from the anti-crossing of the spin and 
orbital states, which are coupled by terms with the x-like 
symmetry. Thus, the anti-crossing will not occur if 



1^]^ cos(7 — (5) — Ijj^ sin(7 + 5) = 0. 



(31) 



Hi=- fiB\\az{x[lgl, cos(7 ~ 5) - l^^ sin(7 + S)] 
+ yi^BR sin(7 ~6)~ l^^ cos(7 + S)]}. 



(30) 



The angles 7 and 6 that satisfy the above equation de- 
fine an easy passage. For a double dot oriented along 
[100] direction {S = 0) the easy passage occurs for an 
in-plane magnetic field oriented along angle 7 given by 
tan 7 = Id/Ibr- Similarly to the single dot case, the 
measured angular dependence recovers the ratio of the 
spin-orbit couplings. Now also revealing which one is 
larger. More important, as can be seen from Eq. H31|l . 
both linear Bychkov-Rashba and Dresselhaus (also cubic) 
spin-orbit terms contribute to the anti-crossing; in single 
dots it is only the Bychkov-Rashba coupling which gives 
relevant spin hot spots. The position of the easy passage 
is then given by an interplay of all the spin-orbit terms. 
If the double dot is oriented along [110] [5 — it /A), the 
condition for the easy passage is 7 = 135°, being inde- 
pendent on the spin-orbit couplings. The importance of 
this result has been pointed already in Ref. |33, where 
the corresponding numerical results are presented. 



B. Perpendicular magnetic field 

1. Orbital relaxation rate 

There are two different regimes for the orbital relax- 
ation, depending on the energy difference of the ground 
and orbital states, E = 26 Et, which is more sensitive to 
the interdot distance than to the confinement length. If 
V = d/ls ^ 1, then E « fi^ /mlg, decreasing with in- 
creasing the magnetic field or the interdot distance. The 
limit of high £ applies and the rates are comparable to the 
single dot case. On the other hand, if I? ^ 1 the energy, 
and thus also the rates, drop exponentially with increas- 
ing the magnetic field or the interdot distance. Due to 
the complex interplay of the magnetic field and interdot 
distance, no power law dependence of the rates on mag- 
netic field can be identified. However, approximations 
in Tab. ^give analytical formulas with a fair agreement 
with numerics, if the energy difference E « 2SEt which 
is given by Eq. (|25J). 

The dependence of the orbital relaxation rate on the 
magnetic field and the interdot distance, for a confining 
length 32 nm, is shown in Fig. [3 The lower left cor- 
ner is the regime of the high £ limit. The rate here is 
similar to the single dot case. The opposite corner is the 
regime of an exponentially small energy difference and 
the rate is practically zero. The transition between these 
two regimes comes for a smaller interdot distance if the 
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FIG. 5: (Color online) Orbital relaxation rate in a double 
quantum dot as a function of in-plane magnetic field for 7 = 
0° and the interdot distance d / tunneling energy SEt, for a 
confinement length 32 nm. The relaxation rate is given on 
the logarithmic scale in the units of s~^. The double dot is 
oriented along [100] (5 = 0°). 



magnetic field is higher, since the transition occurs when 
d ~ Zfl. Again, as in the single dot case, the anti-crossing 
does not have a large influence on the orbital rate - in 
the figure it can hardly be seen. For interdot distances 
much larger than Ib the dots are efl'ectively isolated. 



2. Spin relaxation rate 

Spin relaxation in double dots reveals a surprising com- 
plexity as compared to the single dot case. The com- 
plexity is due to the strong anisotropy of spin hot spots. 
While anisotropy appears already in single dots, caused 
by the interference of the Bychkov-Rashba and Dressel- 
haus couplings, additional anisotropy appears in spin hot 
spots. This anisotropy does not require the presence 
of both couplings. Instead, it is caused by the selec- 
tion rules for spin-orbit virtual transitions in the double- 
dot spectrum. The corresponding physics is described 
by the transformed Hamiltonian Hi of Eq. 1301 We 
have presented the corresponding numerical calculation 
in Ref. |33. Here we discuss the individual contribu- 
tions of the Bychkov-Rashba and Dresselhaus terms in 
the spin relaxation rate and, specifically, in the spin hot 
spot anisotropy. 

The contribution to the spin relaxation rate from the 
Bychkov-Rashba (Dresselhaus) term is shown in the up- 
per (lower) part of Fig. The changes of the upper 
figure, if the Dresselhaus terms were pres ent, would be 
very small (compare with Fig. 2 in Ref. ISTJ) . For low 
magnetic fields the rate grows with increasing magnetic 
field, as we expect from Tab. Q] However, similarly to the 
in-plane magnetic field case, the spin hot spots (ridges 
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FIG. 6: (Color online) Spin relaxation rate as a function of 
perpendicular magnetic field for 7 = 0° and the interdot dis- 
tance d I tunneling energy 8Et (at zero magnetic field only), 
for a confinement length 32 nm. The relaxation rate is given 
in logarithmic scale in the units of s~^. The double dot is 
oriented along [100] (^ = 0°). The upper figure shows results 
when only the Bychkov-Rashba term is present in the Hamil- 
tonian. In the lower figure, only the Dresselhaus terms are 
present. 



in Fig. El) dominate the rate for most of the param- 
eters' range. The interdot distance strongly influences 
the spin relaxation rate by determining the position of 
anti-crossings. In high magnetic fields, the spin state can 
anti-cross higher orbital states depending on the sym- 
metry of these states. However, the influence of these 
anti-crossings on the rate is limited to a narrow region 
of magnetic fields, since the dots are effectively isolated 
at high fields and the crossing states do not comply with 
the selection rules for spin-orbit couplings of single dot 
states. 

It is interesting to compare the contribution to the spin 
relaxation by the Bychkov-Rashba and the Dresselhaus 
terms. Let us first look at the single dot regime, which 
in Fig. iniis visible at d = 0. The spin hot spot appears 
only for the Bychkov-Rashba coupling, in line with our 
earlier observation*. The Dresselhaus coupling becomes 
effective only in the coupled-dot system in which the sym- 
metry of the lowest orbital states allows the coupling at 
the level crossings. The coupling is again absent at two 
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FIG. 7: (Color online) Spin relaxation rate as a function of 
7 and the tunneling energy for _B|| = 1 T, for [110] growing 
direction. The dot orientation is given hy 5 = n/2. The 
relaxation rate is given in logarithmic scale in the units of 



FIG. 8: (Color online) Spin relaxation rate as a function of ^ 
and tunneling energy for _B = 1 T, for [110] growing direction. 
The dot orientation is given by 5 = n/2. The relaxation rate 
is given in logarithmic scale in the units of s~^. 



C. Other growing directions 



isolated dots {d -^ 00). Another nice feature seen in Fig. 
Elis the transformation of the single-dot spin hot spot at 
about 5 T to a double-dot spin hot spot at lower fields, 
while the single-dot spin hot spot that starts at about 9 
T shifts towards 5 T in the double dot and remains there 
at all couplings. 

Similarly to the in-plane field case, we can understand 
the anisotropy of the spin relaxation in perpendicular 
magnetic field by transforming the Hamiltonian of Eq. 
(|22|l into a coordinate system with the x-axes being along 
d: 



Hi ^iiBj^{x[<7x{lg]^ — Ij^ sin 2(5) — (Jylj^^ cos 2(5] 
+ ui'^viJ'BR + ^D^ sin 25) - axlo"" cos 2(5]}. 



(32) 



Thus far we have considered lateral quantum dots de- 
fined in a (001) plane of a GaAs heterostructure. A 
different growing direction leads to a different form of 
the Dresselhaus spin-orbit interactions!: (the form of the 
Bychkov-Rashba term remains unchanged) and to dif- 
ferent conditions for the easy passage. Our results are 
summarized in Tab. IIIII For the [111] growth direction 
the Dresselhaus term has the same form as the Bychkov- 
Rashba one. Our results easily translate for this case by 
placing formally Ijj —> 00. There will be no spin relax- 
ation anisotropy in single dots, while in double dots spin 
hot spots vanish for cos(7 — (5) = at in-plane fields. For 
a general magnetic field a weak passage occurs only at 
specific spin-orbit parameters, given by 2y/3lBB. + lD = 
(the couplings can be negative). 

A less trivial situation occurs for the [110] grown quan- 
tum well. The linear Dresselhaus coupling has the form 



Due to the presence of the orbital effects of the perpen- 
dicular magnetic field, the unperturbed states have no 
specific symmetry under inversions along x. As a result 
only in the limit of low magnetic fields {Is ~ ^0), for 
us below 1 T, the term in Eq. H32|) containing x dom- 
inates over the term containing y; in higher fields both 
terms contribute. In this limit the condition for a sup- 
pression of the anti-crossing is Id = Ibr and 5 = 45°. 
This we call a "weak passage", since the anti-crossing, 
while strongly suppressed, is still present. If the condi- 
tion for a weak passage is not fulfilled, the spin relaxation 
rate, as a function of (5, still has a minimum a,t 5 — 45° 
and a maximum at (5 = 135°. However, the ratio between 
the two extremal values is in general of order 1. 



Hn = — 



Aral 



-CfzP'x 



D 



(33) 



Unlike the Bychkov-Rashba coupling, which has eigen- 
spins always in the plane, the [110] Dresselhaus term has 
eigenspins oriented out of the plane. 

The calculated spin relaxation rate for the double dot 
system oriented along 5 = 7r/2 in an in-plane magnetic 
field of _B| I —1 T is shown in Fig. [T] The spin hot spots 
exist for all orientations of the field except at multiples 
of TT. This is confirmed by analytical considerations sum- 
marized in Tab. IIIII The easy passage exists if the dot 
is oriented along the (rotated) x axis, while the in-plane 
magnetic field is along y. Also, the [110] Hamiltonian 
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growing dir. 


in-plane 


general 


[001] 


Ibr cos(7 + (5) = 
= Id sin(7 — S) 


Id = ^Bfl, (5 = 7r/4 


[111] 


cos(7 — 5) = 


2V3;sH + Zd = 


[110] 


7 = 0, <5 = 7r/2 


^s« cos & = ±2^15 cot ^, 
sin((5 — 7) = ±1 


[cos a sin a 0] 


5 = 7r/2, 
/i3tan7 = — isHC0s2a 


/d = — /sRC0s2a, 
5 = 7r/4, ^ = 



TABLE III: Easy passage conditions for several growing di- 
rections in an in-plane magnetic field and weak passage condi- 
tions for a magnetic field with a nonzero perpendicular com- 
ponent. The z-axis points in the growth direction. The angle 
between d and the (accordingly rotated) a;-axis is 5, the angle 
between the in-plane part of the magnetic field and the a;-axis 
is 7, while ^ is the angle between the magnetic field and the 
z-axis. 



is not invariant under the in-plane inversion of the coor- 
dinates which is why the period in 7 for the relaxation 
rate is twice as in the case of the [001] growing direction. 
However, the part of the Hamiltonian important for anti- 
crossing is invariant with respect to inversion along y. 
Therefore, the results in Fig. [7| for 7 > tt are equal to 
those at 27r — 7 to a very good approximation. 

In order to demonstrate the difference between easy 
and weak passages, we plot in Fig. |Hlthe calculated spin 
relaxation rate in double dots defined in a (110) plane. 
The dots are oriented along y. From Tab. IIIII one gets 
the conditions for the weak passage to be 7 = 0, and, 
for our spin-orbit couplings, ^ = 56°, where ^ is the an- 
gle between the magnetic field and z. This arrangement 
corresponds to the "neck" on the spin hot spot in Fig. 
|H1 However, contrary to an easy passage, here the width 
of the anti-crossing region is finite and gets larger with 
increasing magnetic field (not shown). Since all weak 
passages we found depend on spin-orbit couplings, they 
(better, the corresponding geometries) are much less use- 
ful for robust inhibiting of spin relaxation than easy pas- 
sages. 

In the above analysis we have not considered the cu- 
bic Dresselhaus term, Hdj,, in deriving the conditions 
for easy passages. Being cubic, even after rotating the 
double dot (5 7^ 0), it always has qualitatively the same 
symmetry properties with respect to inversions about x 
and y - it is a sum of two terms, one with symmetry of 
X and one y. Therefore the presence of Hd?, does not 
destroy the easy passage. It can only slightly change 
the conditions for the easy passage to occur. For our 
parameters this change, checked numerically, is only on 
the order of 1° for the of angles in Tab. IIIII so the linear 
terms should provide a realistic guidance to experimental 
demonstrations of the predicted anisotropy. 



V. CONCLUSIONS 

We have calculated phonon-induced orbital and spin 
relaxation rates of single electron states in single and 
double quantum dots. The rates were calculated as a 
function of in-plane and perpendicular magnetic fields, as 
well as a function of the field and (in the case of double 
dots) dots' orientation. Realistic, GaAs defined, electron- 
phonon piezoelectric and deformation potential Hamilto- 
nians were considered. Similarly, relevant spin-orbit in- 
teractions, namely the Bychkov-Rashba and linear and 
cubic Dresselhaus couplings, were used to calculate the 
spin relaxation rate. We have supported our numeri- 
cal findings by analytical models based on perturbation 
theory, deriving effective Hamiltonians which display, in 
the lowest order, all the important effects seen in nu- 
merics. We have proposed using a classifying dimension- 
less parameter E which allows to obtain relevant trends 
and order-of-magnitude estimates in important limiting 
cases. 

In the case of single dots, we have carefully analyzed 
the theoretically predicted anisotropy of the spin relax- 
ation rate in an in-plane magnetic field. The anisotropy 
comes from the interplay of the linear Bychkov-Rashba 
and Dresselhaus terms (if only one of the terms domi- 
nates, the anisotropy is absent). Experimental verifica- 
tion of the anisotropy would give a strong evidence of 
the spin-orbit mechanism of spin relaxation. Further- 
more, such a measurement would enable to estimate the 
ratio of the two relevant spin-orbit terms. 

For single dots in a perpendicular magnetic field, which 
causes cyclotron effects as well as Zeeman splitting, we 
have numerically investigated the orbital relaxation rate. 
In addition, we have provided a simple analytical scheme 
to estimate the rates in the important limits of low and 
high magnetic fields, and found the corresponding rate 
as a function of the confining length. The orbital relax- 
ation rate is found to be of the order of 10^ s^"'^, with 
a relatively small dependence on the magnetic field. At 
anti-crossings the orbital relaxation rate is reduced by a 
factor of two. At low magnetic fields the rate is domi- 
nated by the deformation potential electron-phonon in- 
teraction, while at high fields it is dominated by piezo- 
electric phonons. 

On the other hand, the spin relaxation in single dots is 
always dominated by piezoelectric transversal phonons. 
The contribution of deformation potential phonons is 
more than a decade smaller. The rate is on the order of 
10^ s~^ over a large region of parameters (magnetic field 
and excitation energy). However, the rate is strongly en- 
hanced in the region of anti-crossing/spin hot spot, where 
it becomes comparable to the orbital relaxation rate. We 
have also provided analytical estimates of the rate (away 
from the spin hot spots) for various phonon contribu- 
tions, at the limits of low and high magnetic fields. 

The physics is more complex in coupled dots. We 
have numerically studied spin relaxation in double dots 
in in-plane magnetic fields, in which the rate is strongly 
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FIG. 9; (Color online) Spin relaxation rate in a double dot 
as a function of the orientation of the in-plane magnetic field 
and tunneling energy for B = 5 T, for [001] growing direction. 
The dot orientation is given hy 5 = 7r/4. A small asymmetric 
term is added into the confinement potential (an electric field 
of 10^ V/m in y direction is applied in on one of the dots). By 
this, the easy passa ge i s turned into a weak passage - compare 
with Fig. 4 in Ref.lSTJ 



anisotropic in the direction of both the magnetic field 
and the dots' axis. Similarly to the single dot case, the 
piezoelectric phonons dominate spin relaxation here. We 
have demonstrated that a spin-hot spot exists at useful 
magnetic fields (say, 1 T) and interdot couplings (0.1- 
0.01 meV). In fact, a spin hot spot is a typical phe- 
nomenon in symmetric double dots since it appears when 
the tunneling (coupling) energy becomes comparable to 
the Zeeman splitting. Fortunately, the spin hot spots are 
strongly anisotropic, due to the symmetry of the low- 
est orbital electronic states, and they vanish at certain 
orientations of the field and the dots' axis. We have sys- 
tematically investigated these "easy passages" using an 
analytical model. We have found the criteria for the ab- 
sence of spin hot spots for different growth directions 
of the underlying quantum well. These criteria should 
be seriously considered in fabricating double dot systems 
for spin-based quantum information processing which re- 
quires low spin relaxation. 

For double dots in a perpendicular magnetic field, the 



orbital relaxation rate is most influenced by the energy 
difference of the corresponding coupled states. The en- 
ergy has a range over eight orders of magnitude due to 
the cyclotron effects in the interdot coupling. As in the 
single dot case, both deformation potential and piezoelec- 
tric phonons can dominate the orbital relaxation. The 
spin relaxation in double dots in a perpendicular field 
has similar qualitative features as in the single dot case, 
with an additional anisotropy given by the orientation 
of the double dot with respect to the crystallographic 
axes. However, unlike in in-plane fields, only weak easy 
passages (in which spin hot spots form a neck on the pa- 
rameter map, rather than disappear altogether) exist in 
a perpendicular magnetic field. We have also observed a 
nice shift of spin hot spots to the lower field neighbors 
as the tunneling between the dots decreases. While the 
perpendicular fields provide a nice opportunity to study 
fundamental physics of double dot systems, they are less 
useful in quantum information processing due to the om- 
nipresence of spin hot spots and weak passages. 

Our final note concerns the symmetry of the double dot 
systems investigated in this paper. Do our conclusions 
hold if the symmetry is broken? The answer is yes, if the 
double-dot system still possesses either x- or y-like sym- 
metry. Suppose, for example, that a weak electric field 
is applied along x or y, or one of the dots is somewhat 
smaller than the other. The spin hot spot anisotropy still 
leads to easy passages in spin relaxation in in-plane mag- 
netic fields. On the other hand, if the symmetry breaking 
is xy-like (an electric field pointing along a diagonal, for 
example), the easy passage is destroyed since the selec- 
tion rules for the lowest orbital states will allow coupling 
of the states by the term containing y in Hi of Eq. QUI 
(recall that it was the vanishing of the term containing x 
that lead to the appearance of easy passages). This situ- 
ation is demonstrated in Fig. ^ A double dot system in 
an in-plane field of 5 T is oriented along [110] (the growth 
direction is [001]). If the double dot is symmetric, an easy 
passage exists for 7 = 135° (the corresponding figure is 
given in Ref. [37J). However, if one of the dots is subject to 
a y-like electric field, so that the overall symmetry of the 
perturbation is xy-like, the easy passage turns to a weak 
passage — at all directions of the in-plane magnetic field 
there exists an interdot coupling in which the spin relax- 
ation rate is greatly enhanced. This is another important 
message for spin-based quantum information processing 
in quantum dots. 
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